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Abstract. Total cross sections for electron capture are calculated for collisions of 
fast protons and a-particles with atomic hydrogen. The distorted-wave impulse 
approximation is applied over the energy range 10-1500 keV/u. State-selective results 
are given for the Is, 2s and 2p levels. Both the post and prior forms of the 
model are calculated and compared with results from other theories and experimental 
measurements. In general the model performs very well in comparison with experiment 
over this energy range though discrepancies arise at lower energies. 



1. Introduction 

The study of the interaction between atomic hydrogen and fast ions has an important 
role in the physics of fusion plasmas (Gilbody 1995). It is significant also in other areas 
of research such as the study of stellar atmospheres (Janev et al 1987) and the physics of 
proton aurorae (Van Zyl 1993). An essential requirement for the study of these complex 
phenomena is the availability of accurate data for a wide variety of fundamental ion- 
atom collision processes. Given that the velocity range is such that we can consider 
the nuclear motion as classical and solve the time- dependent equation of motion for the 
electrons, one might suppose that simulations of these processes would be rather simple. 
However this is not the case, and the theory behind this problem remains an area of 
active investigation (Bransden and McDowell 1993). While different strategies have been 
applied to problems of this nature, the most common approach relies on basis function 
expansions (Fritsch and Lin 1991, Kuang and Lin 1996, McLaughlin et al 1997). This 
and other methods, including the the direct solution of the partial-differential equation, 
have been reviewed comprehensively by Bransden and McDowell (1993). 

While the basis function method has a well-founded popularity, it does have some 
drawbacks. Firstly, it can be computationally expensive when excited-state coupling 
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is involved. Moreover the choice of basis set is problematic; one cannot always 
rely on the use of larger basis sets to produce better data. Kuang and Lin (1996) 
indicate some of the spurious effects that can arise if the size and type of a basis 
set are not appropriate. In spite of these caveats, it is fair to say that recent large- 
scale calculations of this type have given excellent agreement with experimental data 
over the appropriate energy range (McLaughlin et al 1997). There exist limitations 
to the method because it is primarily designed for strong diabatic coupling between 
bound states which are nondegenerate. It tends to work less well for higher energies 
at which ionization (continuum coupling) is prevalent, and when mixing with nearly 
degenerate states takes place. While significant advances have been made in the use 
of pseudostates to simulate such couphng (Shm and Ermolaev 1994, McLaughhn et al 
1997) the implementation of the method requires access to high-performance computers. 
An innovation marrying bound and continuum states has been developed by Brown 
and Crothers (1994, 1996a, 1996b) and includes continuum coupling explicitly in the 
wavefunctions. This is known as the symmetrized variational continuum distorted- wave 
(SVCDW) method and it gives excellent agreement with experimental data for electron 
capture, excitation and ionization over the intermediate energy region. 

Our paper is concerned with fundamental charge exchange processes, namely those 
occurring between hydrogen atoms and fast protons or a-particles. Charge exchange 
poses theoretical difficulties in that one requires an accurate description of the electron 
dynamics around both the target and projectile. Often a large two-centre basis function 
expansion is the prescription for achieving this aim, or a restricted optimized expansion 
(Brown and Crothers 1994). Considering alternatives to these approaches, the fact that 
charge exchange rapidly decreases in importance as the energy is increased means that 
one can often use perturbation theory at high energies to obtain reliable predictions 
(Belkic et al 1979). This also has the advantage in being easy to calculate even 
with modest computational resources. However, refinements of perturbation theory, 
including second-Born type corrections, are much more difficult to compute and often 
do not improve the agreement with experiment (Bransden and McDowell 1993). This 
paper investigates the effectiveness of extending one such high-energy approximation 
into an intermediate energy regime. Here we understand the term intermediate energy 
to refer to the range 10-200 keV/u, and this is the range of greatest interest for beam 
injection techniques used in fusion science (Gilbody 1995). The method we adopt, the 
distorted- wave impulse approximation (dwia), though based upon perturbation theory 
can account for high-order potential scattering from both centres. However the theory 
is not of the close-coupling type and thus in its present form cannot describe strong 
intrashell (Stark) mixing or backcoupling (McCann 1992). 

The DWIA was formulated (Miraglia 1982) as an improvement upon the conventional 
plane- wave impulse approximation (PWIA) for electron capture in ion- atom collisions. 
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The PWIA is designed to work for collisions in which the collision is fast in terms of the 
atomic time-scales (projectile or target). The validity for the impulse approximation 
depends upon the collision time being much shorter than some characteristic orbital 
period of the electron bound state. This prerequisite comes from treating the weaker 
interacting nucleus as a spectator particle. If the collision is sufficiently fast on the 
atomic time-scale of this particle then the spectator particle cannot interact strongly 
with the electron during the transition. Its role is confined to describing the electronic 
state before or after the supposed three-body interaction in terms of binding the electron 
and providing it with a momentum distribution in the bound state. This establishes a 
loose validity criterion for the approximation that the collision time should be shorter 
than the internal interaction time, or in other words, the atomic orbital period. In 
applications the impulse approximation was found to be very effective in reproducing 
a large variety of experimental data for very fast collisions (Jakubassa- Amundsen and 
Amundsen 1980). However the inherent asymmetry of the method leads to post and 
prior forms of the transition amplitude; in general these do not agree well. The 
two versions of the approximation converge at high energies, but strongly diverge at 
intermediate and low energies. On the other hand the DWIA seems to be better suited to 
both symmetric and asymmetric collisions. While we find that there is still a post-prior 
discrepancy, it is very much less than that which arises from using the PWIA. We employ 
this model to investigate the following coUision processes 



These are ideal processes to test any theory because the atomic structure is trivial. Also 
the absence of screening effects allows us to write continuum states in terms of Coulomb 
functions. For this reason the evaluation of matrix elements is straightforward. 

2. Method 

In this three-body re-arrangement reaction, we have a bare projectile nucleus of charge 
Zp travelling at a velocity v relative to the target atom which has nuclear charge Zt- 
The position of the electron in question, with respect to the target and projectile nuclei, 
is defined to be and Vp, respectively. The initial and final atomic energies of the 
electron are denoted by Et and Ep, corresponding to the target and projectile state 
electron wavefunctions and $/. The momentum-space representations of these states 
are labelled as and $/ in this paper. 

The definitions of post and prior can lead to confusion so let us clarify their meanings. 
We use the standard convention that post refers to using the impulse approximation in 



H+ + H(ls) ^ H(ls, 2s, 2p) + H 



(1) 



and 




(2) 



the initial channel for the interaction between the projectile ion and target atom. In 
order to calculate the prior form, a separate program must be written. 

DWIA(post) P(PWIA) + (e-, T) ^ (P, e") + T(cDW) (3) 
DWIA(prior) P(cDW) + (e", T) (P, e") + T(pwiA) (4) 
In the prior form, the DWIA transition amplitude a^j^\b) for a given impact parameter 
b consists of a continuum distorted-wave ^ctdvf entrance channel and a distorted- 

wave x^iA from the PWIA in the exit channel: 

a\-\b) = -^7dt {xU I H - id, I ^Sh.) (5) 

— oo 

The expression in momentum-space (McCann 1992) in the prior form is: 

S^T)(ry) = {'l'n)-^l'^v-^iN{ap) J dk N{aT) ^{k -v)c d (6) 

c = j drr exp(i/x • r^) iPi [iar, 1; i(/crr + • r^)] ^/^^'^'^(^^t) (7) 
d = J drp exp{iu: • rp) Vrp iFi[iap;l;i{vrp + V • rp)] (8) 

with the momentum transfer vector /x = —77 + {AE/v"^ — ay = Zx/k, ap = Zp/v 
and oj = —fj, — k, the electron momentum being k. The transverse impulse is denoted 
by rj. The normalization factor N{x) is defined by N{x) = exp(7ra;/2) r(l — ix), and 
AE = Bp — Erp- The post form of the transition amplitude is: 

~a^ff{r]) = {27T)-^/\-HN{a'T) J dk N{a'p) ^^{k + v) c' ■ d' (9) 
c' = j drp exp{i^i' ■ Vp) iFi [ia'p] 1; i{krp -k-Vp)] Vrp^){rp) (10) 

= y drp exp(icj' • r^) Vj.^ iPi[ia^; 1; i{vrT + f • t^-)] (11) 

Here the momentum transfer vector fj,' = rj — [AE/v"^ + ^)v, a'p = Zp/v, a'p = Zp/k 
and (jj' = —fx' + k. 

For the purposes of our calculation we have considered transitions between excited 
target states and the projectile ground state. One can always relate the data obtained to 
the inverse reaction using the principle of detailed balance and time-reversal symmetry, 
i.e. cross sections for the processes ([^) and (0) can be elicited by looking at their 
time-reversed counterparts: 

H+ + H(ls,2s,2p) ^ H(ls) + H+ (12) 

and 

H+ + He+(ls, 2s, 2p) ^ He^^ + H(ls). (13) 

In our calculations we have therefore explicitly considered electron capture into the Is 
state from the Is, 2s and 2px,y^z states. From now on, we will discuss calculations in 
terms of the prior form only, for the sake of convenience. 



5 



2.1. 2s-ls transition 

Previous work using the DWIA (Gravielle and Miraglia 1988) dealt with electron capture 
from the K shell, specifically Is-ls electron transfer. It was noted (McCann 1992) that 
the evaluation of the azimuthal integral, with the polar axis lying along the momentum 
transfer vector /x, generates an Appell function Fi (Appell and Kampe de Feriet 1926) 
which is a hypergeometric function of two variables. Its many linear transformations and 
analytic continuations (Olsson 1964) can be used to evaluate this function. Nevertheless 
it remains a major obstacle to efficient and accurate computation of the capture cross 
section. 

For the case of 2s-ls electron transfer, it can easily be shown using parametric 
differentiation that the spatial integrals, (|^) and (§), over vt and rp simplify in a manner 
very similar to that found in the Is-ls symmetric transfer calculation (McCann 1992). 
Using spherical coordinates for the fc-integral in (^, we denote the angles between ^ 
and k, and that between /x and v by 6^/^ and 6^^. Taking the angle between these two 
planes as 0, we then encounter the following term: 

J d0 {A + B COS (Py^C + D COS (P)-^-'^'' = 27T {A + B)-\C + D)-^-'^'' 



xFi(|; 2,1 + zap; l;xi, 0:2) (14) 

with xi = 2B/{A + B) and X2 = 2D/{C + D). The quantities A, B, C and D have the 
definitions: 

A = Zp + v'^ + k"^ — 2vk cos 9^k cos 9^^ 

B = 2vk sin 9^k sin 9^^ 

C = —AE + v'^ /2 — vk cos 9^^ cos 9^^ 

D = vk sin 9^}; sin (15) 



2.2. 2p-ls transition 

The lack of spherical symmetry for the 2p sub-states means that the choice of the axis of 
quantization is important (Coleman and Trelease 1968). Using the momentum transfer 
vector fi as the polar axis simplifies calculations and the change to the laboratory frame 
of reference requires a simple rotation transformation. 

A different line of approach takes advantage of parametric differentiation. We can 
consider the spherical harmonics in their real form: 

^2p,irT) = i4/'(27r)-i/2exp(-ATrT) ttj (16) 
where At = and j G {x, y, z}. In evaluating (|^, integration by parts gives: 

C = /l+/2 (17) 
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Table 1. Cross sections ani (10~^^cm^) of electron capture for proton-hydrogen 
collisions at energy E keV: Is to nl transition. The integer in parenthesis indicates the 
power of ten by which the number has to be multiplied. 



E(keV) 




125 


250 


500 


750 


1000 


1500 




DWIA 
CDW 


2.45(-l) 
2.70(-l) 


1.35(-2) 
1.38(-2) 


4.83(-4) 
4.83(-4) 


5.97(-5) 
5.96(-5) 


1.30(-5) 
1.29(-5) 


1.43(-6) 
1.43(-6) 


0'2s 


DWIApo.t 
CDW 


4.68(-2) 
4.50(-2) 
5.05(-2) 


2.19(-3) 
2.20(-3) 
2.26(-3) 


6.98(-5) 
7.02(-5) 
7.13(-5) 


8.23(-6) 
8.28(-6) 
8.43(-6) 


1.74(-6) 
1.75(-6) 
1.79(-6) 


1.87(-7) 
1.88(-7) 
1.93(-7) 


CT2p 


DWIAp,,,* 
CDW 


9.89(-3) 
1.29(-2) 
1.54(-2) 


2.65(-4) 
3.76(-4) 
4.83(-4) 


4.93(-6) 
7.38(-6) 
1.01(-5) 


4.42(-7) 
6.60(-7) 
9.53(-7) 


7.78(-8) 
1.15(-7) 
1.74(-7) 


6.54(-9) 
9.67(-9) 
1.58(-8) 



-J dvT exp{in • Tj) $j(rT) V^^, iFiliar, l;i{krT + k • Vt)] (18) 
-ifj, J dry exp(i/x • rj.) iFi[iaT; 1; i{krT + k • r^)] 



(19) 



The integral /i is further simphfied by writing it as: 
d 

I drT — ^— ^-^ —<PAr-r]il'i hai-; 1; i[qrT + q • T^-jJ 

q=k 

exp(— At?"t + i^J' • fx) 



Ii = -q 



dq 



d d 



dq dfij 

XiFiliar, 1; i{qrT + q • Vt)] 
d d 



drn 



TT 



1-1/2 



iq 



dq dfij 



q=k 



laT 



q=k 



(20) 



where a = |(Ar^ + A*^), (3 = fj, • q — iXx q, and d/dfij calls for partial differentiation 
with respect to the j component of /j,, j G {x,y,z}. The momentum transfer vector 
/X consists of two mutually orthogonal components, with one running parallel to the 
direction of travel v of the incoming projectile which we can take as the 2;-axis. The 
transverse component is —r] and this lies within the xy-plane. It is then convenient to 
let r] lie along the x{y)-a.xis with the result that capture from the 2py(^x) state is parity 
forbidden. 



Employing the same method, the expression (19) can be shown to be given by: 



lV/2(2vr)-i/V 



d d 



27ia'^^-^ 
{a + f^Y'^T 



(21) 
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For electron capture from the 2p states, the azimuthal integral of (0) has the form: 



2lT 



J = Jd(j){A' + B' cos (j))-^ [C + D' cos 0)-^-^'^^ {E' + F' cos 0) (22) 



Like those of (|T3p, the quantities A' , B', C, D', E' and F' depend upon r], k and the 
polar angle between k and /x, but are much more complicated; we do not quote the 
details. 

Parametric differentiation transforms (p2D into an expression containing derivatives 
of the more familiar Appell functions Fi. 



Since 



^ OA' dB''^ V ; V -r ; 

X Fi(i;l,l + mp;l;xi,X2)} (23) 



d^^Fi{a;hi,h2;c;xi,X2) = —Fi{a + l;bi + l,b2;c+ l;xi,X2) (24) 

c 



we therefore have: 

J = ~2n{A' + B'y^C' + D')"'"'"^{i^i(i; 2, 1 + zap; 2; Xi, X2) (A'F' - 5'E') 

-Fi(i;l,l + mp;l;xi,X2) (A' + 5')(^' + F')} (25) 

As with Is-ls electron capture, the corresponding 2s-ls and 2p-ls transition amplitudes 
possess two singularities, at C±D = and C'±D' = respectively, which are integrable. 

2.3. Symmetric models 

In the DWIA we inevitably have a post-prior asymmetry. If there is an inherent 
asymmetry in the reaction, for example Zt ^ Zp, the choice between post and prior 
models can be argued in favour of one or the other. A detailed discussion of this 
question is given in the next section. A merit of the CDW theory is the intrinsic post- 
prior equivalence. Symmetry can be artificially introduced in our model by averaging 
the post and prior amplitudes. The resulting symmetrized DWiA is given by: 



4(^) = l[4^(^) + 4tH^)J (26) 

Miraglia (1982) had also proposed making the impulse approximation in both entry and 
exit channels (generalized impulse approximation gia). While it is much more difficult 
to calculate the GIA, it is not clear whether this GIA model will have a larger range of 
validity than the DWiA, and thus whether it offers a substantial gain in accuracy. 
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Table 2. Cross sections dni (lO^-'^^cni^) of electron capture for helium nucleus- 
hydrogen collisions at energy E keV: Is to nl transition. The integer in parenthesis 
indicates the power of ten by which the number has to be multiplied. 



E(keV) 




10 


25 


60 


100 


150 


200 




BWlApost 


2.76(0) 


2.40(0) 


9.02(-l) 


4.87(-l) 


2.42(-l) 


1.25(-1) 






3.66(0) 


4.22(0) 


2.10(0) 


8.58(-l) 


3.24(-l) 


1.41(-1) 




CDW 


2.78(1) 


3.07(0) 


2.73(0) 


1.10(0) 


3.94(-l) 


1.64(-1) 




SVCDW 


2.92(-l) 


9.93(-l) 


1.70(0) 


7.67(-l) 


2.86(-l) 


1.21(-1) 




BWlApost 


2.07(1) 


4.81(0) 


1.30(0) 


4.62(-l) 


2.54(-l) 


5.42(-2) 




DWIAp„o^ 


1.77(1) 


5.86(0) 


1.33(0) 


4.47(-l) 


1.38(-1) 


5.11(-2) 




CDW 


1.68(2) 


2.78(1) 


3.53(0) 


7.84(-l) 


1.94(-1) 


6.47(-2) 




SVCDW 


1.58(1) 


9.41(0) 


1.90(0) 


4.72(-l) 


1.25(-1) 


4.44(-2) 




BWlApost 


1.05(2) 


1.46(1) 


2.26(0) 


6.21(-1) 


1.53(-1) 


4.52(-2) 






1.32(2) 


3.54(1) 


5.79(0) 


1.02(0) 


1.93(-1) 


5.20(-2) 




CDW 


1.00(3) 


1.10(2) 


8.01(0) 


1.23(0) 


2.24(-l) 


6.00(-2) 




SVCDW 


6.07(1) 


2.46(1) 


2.95(0) 


5.32(-l) 


1.12(-1) 


3.18(-2) 



3. Results and discussion 

We present results for capture cross sections as a function of energy and compare this 
data with experiment and other models. For the most part we use logarithmic graphs, 
but in addition we have compiled a small sample of the data in tabular form for reference 
purposes (tables 1 and 2). This helps to form detailed comparison of the post and prior 
forms of the DWIA results which would not be apparent from the graphs. We have 
used the CDW model as a benchmark (Belkic et al 1979) for the results. This model 
reproduces electron capture data very well over the high-energy range but suffers from 
a lack of unitarity at lower energies and invariably produces gross overestimates for 
capture at intermediate energies. 

In figure 1 the DWiA results for electron capture cross sections of the ls-2s transition 
in reaction (|l|) are shown. The corresponding data of Brown and Crothers (1996a) have 
not been presented, although they follow the experimental results very closely, even 
down to 10 keV/u. Instead our results are compared with experimental data and the 
continuum distorted-wave (CDw) results over the intermediate energy range: 20-250 
keV/u. At the upper end of the energy regime we note that both DWiA curves merge 
with the CDW results. However this convergence is not uniform. If we refer to table 1 
for energies beyond the scope of the graph, apart from noting the sharp fall in the size 
of the cross section (o"2s), it is clear that at very high energies the post and prior results 
converge but begin to depart from the CDW data. This can be understood from noting 
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the importance of second-order effects at tliese energies. The CDW theory can partly 
but not fully account for these effects, and the DWIA data are more reliable estimates 
in this case. Nonetheless the CDW results are more than adequate for good estimates of 
cross sections at high energies. 

The intermediate energy range covered by figure 1 indicates the disparity of these 
models more clearly. The CDW results continue rising sharply as is well-known (Belkic et 
al 1979). This is typical of the difficulties in applying high-energy perturbation theory 
over the intermediate energy range. We note that DWiA gives a slight improvement over 
CDW below around 50 keV/u, although its accord with experiment remains poor at lower 
energies. Furthermore the post and prior DWiA models show large differences below 50 
keV/u. In common with the PWIA, the inherent post-prior asymmetry is expected to be 
greatest at the lowest velocities. As a general rule, the discrepancy is worse for PWIA 
than for DWIA (Ng and McCann 1997). 

On theoretical grounds we can propose a prescription for the energy range of validity 
of the DWIA and for the preferred form {post or prior) of the theory in this particular 
reaction. Using the peaking criterion (McCann 1992), the orbital velocity of the electron 
in the final state is smaller than that in the initial state, and so one should apply CDW 
refinements in the exit channel and allow the PWIA to take account of momentum spread 
in the entrance channel. This implies that the post form should be preferred. Next, one 
can use the impulse hypothesis (Gravielle and Miraglia 1988) to establish the energy 
limit - given the largest orbital velocity is 1 a.u. the desired stipulation would be that 
E ^ 25 keV/u. This inequality is consistent with the experimental data in figure 1. 
Since we wish to know the practical limits of our model, the important question here is: 
how strong is the inequality? In this paper we aim to answer this question by a mixture 
of comparisons with experiment and other theories. 

We have also computed, for the sake of completeness, several results using the 
symmetrized DWiA (equation ^6]). The results are shown in figure 1 and give the best 
agreement of all the calculations. Nonetheless this averaging procedure seems rather 
artificial to us, and the presence of a large post^prior discrepancy at the lower energies 
indicates that neither the post nor the prior form is very satisfactory around 25 keV/u. 
On the whole the theory of Brown and Crothers (1996a) still appears to be the most 
satisfactory description of the process both in physical and quantitative terms. 

Figure 2 shows capture cross sections for ls-2p transition. Below 40 keV/u, DWiA 
again gives slightly better results than CDW. In the high-energy range the results do 
not converge. Both sets of DWIA results differ and they in turn are lower than the CDW 
results by a fairly constant ratio (table 2). Again the post DWiA is expected to be the 
best physical model for this process and this is confirmed by the better agreement with 
the SVCDW of Brown and Crothers (1996a) over the energy range 20-100 keV/u. It has 
been noted that the experimental results all have systematic uncertainties with regard 
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to the normalization of the data. However these uncertainties are not large and a series 
of experiments has been carried out on this reaction process using independent estimates 
of the absolute magnitudes of the cross sections. It is very likely the experimental data 
in the figure are reliable both in terms of the energy dependence and absolute values. 
In any case we have not attempted to renormalize the data to theoretical predictions. 

The process (|^) is dominated by Is and 2s transitions above 100 keV/u. For 
higher energies we can estimate capture into higher excited states (cr„) by the rough 
approximation based on the n~^-distribution of populations. The total cross section 
therefore becomes atotai = o"i + 1.62(72 where ui and a2 are cross sections of capture into 
the n = 1 and n = 2 levels respectively. Results for DWIA and CDW, in the energy range 
125-2500 keV/u (figure 3), are very similar and are in good agreement with experimental 
data. 

The proton-hydrogen collision calculations do not contain much post-prior 
discrepancies due to the nature of the charge exchange process being more or less 
symmetric. Their agreement with experimental data is very satisfactory. However the 
study of reaction (Q) which has a strong in-built asymmetry would provide a much clearer 
critical assessment of DWiA. But first, the nature and implication of the asymmetry must 
be considered with care. 

Generally in electron transfers from the ground state to one which is excited, an 
important point to consider is the momentum spread of the electron, before and after 
the collision. A sharply peaked momentum spread could be treated more efficiently 
by CDW with its peaking approximation. The impulse approximation should then be 
applied to the electronic state with the less sharply peaked momentum distribution. 
The momentum spread can be affected by the electron's distance from the nucleus and 
also by the magnitude of the nuclear charge, and is roughly of the order of ~ Z/n. 
Since the only asymmetry in proton-hydrogen collision comes from the difference in 
initial and final states and not from projectile and target nuclear charges, its post-prior 
discrepancies are not appreciable. The same, however, cannot be said for (|^). In the 
forward reaction, we would expect an electron attached to the heavier (projectile) a~ 
particle to have quite a broad momentum distribution. Since this situation arises after 
the collision, it would seem reasonable to use the prior form of DWiA. The prior form 
corresponds to a PWIA wavefunction in the exit channel, making the bare proton of 
the hydrogen atom more of a spectator particle than the helium nucleus. It has been 
argued (Gravielle and Miraglia 1991) that the converse should be true, the argument 
being based upon the idea of associating the impulse approximation with the stronger 
potential (i.e. in the entrance channel). We take the opposite viewpoint as the weight of 
evidence presented in this paper lends it sufficient support. Consequently for the time- 
reversed reaction (p^, we apply inversion so that its post form becomes favourable. 

Cross sections for selected energies for reaction (2) are tabulated in table 2, the 
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energy range (7-200 keV/u) being the same as that covered by figures 4 and 5. These 
show cross sections of capture into the 2s and 2p states. Comparison is made with 
the symmetrized variational continuum distorted-wave (SVCDW) method (Brown and 
Crothers 1996b), the conventional CDW theory and experimental results. In the case of 
the 2s transition DWIA fares rather poorly. While the curves follow the SVCDW curve 
at higher energies and avoids the low-energy divergence of the CDW, the agreement 
with experiment is not good over the important plateau region (figure 4). The DWiA 
completely fails to describe this feature. Results for the 2p transition (figure 5) on the 
other hand are less conclusive. The shortcomings of the DWiA around 50 keV/u for the 
2s results could be attributed to the lack of strong couphng between the resonant states 
(H-ls, He'''-2s, 2p). Although the DWiA underestimates the 2s cross sections (figure 
4) and overestimates the 2p results (figure 5), this can be explained by a second-order 
process involving the redistribution of populations: H{ls) He^{2p) — > He~^{2s). The 
inclusion of this effect in the SVCDW model brought good agreement with experiment. 
Its neglect in the DWIA means that there is a disconcerting gap between theory and 
experiment in figure 4, with a fraction of the cross section that should have gone towards 
the 2s transition ending up in the 2p instead. 

The total cross section for capture into all states is shown in figure 6. The DWiA 
prior curve, with the scaling law included, performs remarkably well against both 
SVCDW and CDW. Major contribution to the total cross section comes from capture into 
the 2p state and DWIA models this quite well in its prior form, yielding good results. 
The use of the Is, 2s, 2p results to extrapolate for capture to all states is doubtful, of 
course. A simple test shows that when the scaling law is applied to obtain CDW 
results for (T„ (n > 2), it underestimates the true CDW o"„ in the energy range 7-200 
keV/u. This has the effect that the extrapolation can be in error by as much as 60% . In 
view of this, DWiA results for total cross sections should be considered as underestimates 
for capture to excited states and this could explain part of the shortfall between theory 
and experiment that we see in figure 6. If we were to assume that the underestimation 
for CDW will be of the same order of magnitude for DWiA at a fixed energy, appropriate 
augmentation of the results would bring them quite close to the experimental values. 

4. Conclusions 

In conclusion, despite its inherent post-prior discrepancies, DWIA can perform well even 
at energies as low as 20 keV/u if sufficient care is taken in choosing whether to use 
the post or prior form of the theory. This is very useful in establishing the energy 
range of validity. However the computational effort in evaluating DWiA is much greater 
than either CDW or PWIA model; hours as opposed to seconds on one of our RS/6000 
workstations. The additional effort brings its reward in terms of better estimates of cross 
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sections towards the upper end of the intermediate energy range. This does provide very 
useful data for apphcations. There are awkward features of the calculation that make 
it difficult to compute without careful attention. However the method is significantly 
less expensive to compute than large-scale close- coupling methods, though it does not 
give the same accuracy or reliability over the intermediate energy regime. It is difficult 
to extend the DWIA method further. One might attempt a close-coupling approach akin 
to the SVCDW method. We take the view that significant improvements are difficult 
within the DWiA framework and that instead the SVCDW approach seems to combine 
the correct physical and mathematical features of the charge transfer processes in an 
elegant manner. This approach seems the most promising avenue for future work. 
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Figure Captions 

• Figure 1: Cross sections for Is — >• 2s capture for process (|l]): DWIA(prior) ( ), 

DWIA(post) ( ), symmetrized DWIA (o), CDW ( ) and experimental results 

by Morgan et al (1980) (•). 

• Figure 2: Cross sections for Is — > 2p electron capture for process (|1|): DWIA(prior) 

( ), DWlA(post) ( ), CDW ( ) and experimental results by Barnett 

(1990) (•). 

• Figure 3: Total cross section of electron capture into all states in collisions between 

proton and hydrogen: DWIA ( ), CDW ( ) and experimental results by 

Toburen et al (1968) (•), Gilbody and Ryding (1966) (A) and Schwab et al (1987) 

(□)• 

• Figure 4: Cross sections for Is 2s electron capture for process (0): DWIA(prior) 

( ), DWlA(post) ( ), CDW ( ), SVCDW ( — ) and experimental results 

by Shah and Gilbody (1978) (•). 

• Figure 5: Cross sections for Is 2p capture for process (|l]): DWIA(prior) ( ), 

DWlA(post) ( ), CDW ( ), SVCDW ( — ■ — ) and experimental results by 

Ciric et al (1985) (•). 

• Figure 6: Total cross section of electron capture into all states in Collisions between 

helium nuclei and hydrogen: DWIA(prior) ( ), DWIA(post) ( ), CDW ( ), 

SVCDW ( — ■ — ) and experimental results by Shah and Gilbody (1978) (•). 
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